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Abstract. We obtain large n asymptotics for products of powers of the absolute values of 
. the characteristic polynomials in the Gaussian Unitary Ensemble of n x n matrices. Our 
results can also be interpreted as asymptotics of the determinant of a Hankel matrix whose 
symbol is supported on the real line and possesses power-like (Fisher-Hartwig) singularities. 
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1 Introduction 

In the present paper we compute large n asymptotics for the following averages of the 
characteristic polynomials over the Gaussian Unitary Ensemble (GUE) ofnxn matrices H: 



X 



3=1 

where 



I det(g - )gue = ^^'.'.^f ' Ua i>-1> ^ = 1 "--' m ' C 1 ) 



poo i' eye n " L 

D n (a 1 ,...,a m ) = — •••/ J\( Xi _x i) 2 JJ w(x fc )da; fc , w(x) = JJ \x-fij\ ^ . 

(2) 

and all fij are distinct. The condition dtotj > — | guarantees convergence. 
It is easy to verify that 

/oo 
x k w(x)dx, (3) 
■oo 

that is, D n (oei, . . . , a m ) is a Hankel determinant with symbol w(x). 
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The determinant D n (0, . . . , 0) is an n-dimensional Selberg integral [29], it can be evalu- 
ated exactly, and its large-n asymptotics readily can be written (cf. |36j): 



ra-1 



D n {0, . . • , 0) = (27r)"/ 2 2-" 2 / 2 J] j\ = (27r) n (n/2) n2 / 2 n- 1 / 12 e- (3/4)n2+c ' ( - 1) (l + 0(l/n)), (4) 

3=1 

where ('(x) is the derivative of Riemann's zeta-function. 

The average (JTJ is related to several interesting questions. We mention a conjectured 
relation of (JTJ) for m = 1 and such averages over other ensembles of random matrices to the 
mean values {1/T} f Q |C(l/2 + it)\ 2a dt of the zeta-function on the critical line [23], |2H E]- 
Also note that D n (a\, . . . , a m ) is a Hankel determinant with symbol on R and Fisher-Hartwig 
singularities at x — fij. Although asymptotics of Toeplitz determinants with general Fisher- 
Hartwig symbols have by now largely been determined (at least the leading terms thereof, 
see, e.g., [17]). the Hankel case is still open. 

For connections of (pQ) to the one- dimensional gas of impenetrable bosons and other phys- 
ical problems, see [I7J E3 E3 [20] and below, (see also [3], ©, [20], [30], [31], [33] regarding 
subleading asymptotic terms, negative moments and other random matrix ensembles). 

The asymptotics of averages (CQ) for otj positive integers have been found by Brezin and 
Hikami [7J, Forrester and Frankel [T7J, Garoni [21]. For such a/s, the expression (EJ) can be 
reduced to the Hermite polynomials and their derivatives at the points fij. (In the framework 
of the present paper, this can be shown with the help of Christoffel's formula (e.g.. [32], p. 29), 
cf. pj) It is not, however, the case for noninteger ct,'s. We prove 

Theorem 1 Fix Xj and ctj such that Xj e (—1, 1), Xj ^ Xk, for j ^ k, and §laj > —1/2, 
j,k = 1, . . . , m. Then, as n — > oo ; 



(H i det(# - a^i^^gue = n - x )T ]i2 (2) exp i( 2A ? - 

l + O 



n 

l<i<j<m 

where 



Inn 



n 



C(a)=T(a + l/2r 2a exp^2 J In T(s + l/2)cto + a 2 j = 2 W , (6) 

and G(z) is Barnes ' G-function. In |3J), the values of the roots positive for real ctj ( and their 
analytic extension for complex atj) are taken. 

With increasing effort, one can compute an arbitrary number of terms in these asymp- 
totics (and verify a heuristic calculation of them by Gangardt [20] for m = 1). 

The expression ([5]) without the error term was conjectured in [TTJ, Conjecture 4, as an 
extension off integer values of atj. 



(5) 
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The asymptotics for the Hankel determinant D n (ai, . . . , a^) are now just a combination 
of (PQ)) P|) . and ([5]). The corresponding asymptotics for a Toeplitz determinant with symbol 
(on the unit circle) having such type of singularities (multiplied by a rather general "regular" 
function) were found by Widom [37] . The multiplicative constant in the asymptotics involves 
the same combination of G functions as here. Note that averages of type (CD) for the Circular 
Unitary Ensemble are such singular Toeplitz determinants. The conjectured relation to 
Riemann's zeta-function averages was first formulated by Keating and Snaith [23J for that 
ensemble. 

Let m = 1. Since the density of the scaled eigenvalues in the GUE is given by ip(X) = 
(2/7r)\/l — A 2 , — 1 < A < 1 (the Wigner semicircle law), we have that D n (a) is proportional 
to ip a . Similar dependence on zero density as well as presence of the factor G(a + 1) 2 /G(2a + 
1) is conjectured for the large T asymptotics of the zeta-function average mentioned above. 

For a = 1/2, m = 1, Theorem 1 gives an analogue for the GUE of the average of 
the absolute value of the characteristic polynomial found by Fyodorov for the Gaussian 
Orthogonal Ensemble [19] (in the Unitary Ensemble, his method would correspond to the 
case a = 1). 

The condition A G (—1, 1) means that A is inside the bulk of the scaled spectrum of 
H (in other words, inside the support of the equilibrium measure for the GUE). The case 
A G C \ [—1, 1] corresponds to a Hankel determinant with symbol "regular" on the support 
of the equilibrium measure, and the average ([T]) is then found as a particular case of a result 
by Johansson [22] , although some nonessential assumptions in [22] need to be relaxed. (This 
result gives a Hankel analogue of the Szego asymptotics for Toeplitz determinants with a 
regular symbol. The need to consider the equilibrium measure is an additional difficulty in 
comparison with the Toeplitz case.) For real A, a, m — 1, we obtain from [22] as n — > 00: 

/|\|+ A 2 _ -1 \ 2a ( n+a ) 
(| det(H - Ax/2^)| 2Q ) GUE = (2n)«" (A 2 - I)-* 2 H + v o x 



V 7 (7) 

exp 1 2cm (V- \\\VX^l-^ [l + o(l)], AGR\[-1,1], a>-l/2. 

This (including the error terms) and the asymptotics for the border case of A close to ±1 
could also be computed by our methods. It should also be possible to generalize Theorem 
1 replacing x 2 in w(\/2nx) by a more general function V(x). Thus it should be possible to 
verify Conjecture 8 of Forrester and Frankel \17\. 

The one-body density matrix p(Ai,A2) for the ground state wave function of the one- 
dimensional gas of impenetrable bosons equals [IS] , up to a simple factor, to D n (l/2, 1/2) 
(i.e., the corresponding w(x) = \x — X\\\x — A2 |e~ x ), where n represents the number of 
particles. The eigenvalues of the operator with kernel p(\i, A 2 ) are interpreted as occupation 
numbers of effective single-particle states of the Bose system. To prove the expected asymp- 
totic ^^-proportionality of the leading occupation numbers in this one-dimensional system 
(note that Bose-Einstein condensation, which is not expected in one dimension, would cor- 
respond to the n-proportionality) , we need to analyze asymptotics of _D n (l/2, 1/2) for all 
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values of Ai, A2. Theorem 1 gives the desired asymptotics (conjectured in [IS] , eq.(105)) for 
— 1 < Ai < A 2 < 1, but we still need to analyze, in particular, cases of Ai, A 2 approaching 
each other and the points ±1. This will be the subject of a subsequent publication. 

We prove Theorem 1 using Riemann-Hilbert problem methods (see [8] for an introduc- 
tion). Our approach is as follows. Consider the polynomials p n (x) = >c n x n + • ■ ■ orthonormal 
w.r.t. to the weight w(x) = YlJLi \ x ~~ ^jV2n\ 2aj exp(— x 2 ): 

/oo 
Pj(x)p k (x)w(x)dx = 5 jk , j, k = 0, 1,... (8) 
-00 

Then, 

n-l 

D n (ai, . . . ,a m ) = Y\_ x 7 2 ' ( 9 ) 
3=0 

by a well-known formula [321 [8] (for real <x,-; the case of complex otj will be explained below). 
The asymptotics of the polynomials p n (x) (for A = and m = 1) as n — > 00 were recently 
analyzed by Kuijlaars and Vanlessen in [27] by a Riemann-Hilbert problem approach. Using 
|27j . one can obtain any number of asymptotic terms for p n (x) and x n . 

Thus, knowing Xj, j — > 00, we can shed light on the asymptotics of D n (a>i, . . . ,a m ). 
However, since the product of the first Xj, j = 1, 2, . . . remains unknown we are bound to 
lose at least the constant (in n) factor in such obtained asymptotics of D n (cf. [28]). This 
problem is avoided here by deriving an identity for (d/da u ) ln£) n (ai, . . . , a m ), v = 1, . . . , m. 
(cf. Deift's formula [5] and differential identities of [H]), an idea which has proven useful 
in similar situations [T5J EJ [251 E31 El- It turns out that the particular structure of the 
weight w(x) allows us to write the above logarithmic derivative only in terms of p n (z), 
JT OQ Pn(x)w(x)dx/(x — z), and similar expressions with index n — 1, which are precisely the 
quantities whose asymptotics (uniform for a in a bounded set) follow from the Riemann- 
Hilbert problem for p n (z). A fact that considerably simplifies calculations is that these 
quantities only need to be evaluated at a finite number of points z (fij and infinity); in this 
sense the identity we obtain is "local". (As discussed in Section 3, a "nonlocal" identity 
involving integrals with p n {x), p n -i(x) exists for general weights). The identity is found in 
Section 3. Setting in it all ay = 0, v = 1, and integrating over a% from zero to some ai, we 
obtain the asymptotics for hi(D n (ai, 0, . . . , 0)/D n (0, . . . , 0)). Now fixing in the identity a±, 
setting a 2 = ■ ■ ■ = ot m — 0, v = 2, and integrating over a 2 from zero to some a 2 , we obtain 
ln(D n (a 1 , a 2 , 0, . . . , Q)/D n (a\, 0, . . . , 0)). Continuing this procedure, we prove the theorem 
by induction. Note that instead of zero, we could have used any positive integer a (where, 
as noted above, the asymptotics of D n are known) as a starting point for the integration. 



2 Riemann-Hilbert problem for p n {z) 

Consider the following Riemann-Hilbert problem for a 2 x 2 matrix valued function Y(z) 
and the weight w(x) = YYjLi \ x ~ ^jl 2 ^ exp(— x 2 ), 3?«j > —1/2: 
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(a) Y(z) is analytic for z G C \ R. 



(b) Let x G M \ U^ 1 {yU J }. Y has continuous boundary values Y+(x) as z approaches x 
from above, and Y_(x), from below. They are related by the jump condition 

Y + (x)=Y-(x)(l W [ x) y x£l\U™ =1 W. (10) 

(c) Y(z) has the following asymptotic behavior at infinity: 

y-(*) = (j + oQ))(£ asz-oo. (11) 

(d) The matrix of Y(z) is 0(1) for 9?<x,- > 0, and ( ^Szl W/\ \2a \ J f° r < as 
z — > fMj, j — 1, . . . , m, z G C \ R. 

(Here and below 0(a) stands for 0(|a|).) 

It is easy to verify that, provided the system of orthogonal polynomials Pk{z) = XkZ k + 
• • •, xt 7^ 0, k = 0, 1, . . ., satisfying (JSj) exists, this problem has a solution given by the 
function: 

Y( Z )= n ' n J -°° iTTi \ ,, 2 ) 



For general weights this fact was noticed by Fokas, Its, and Kitaev [16] and, in conjunction 
with the steepest descent method of Deift and Zhou [10] for asymptotic analysis of matrix 
Riemann-Hilbert problems, has allowed in recent years to solve many previously unaccessible 
asymptotic problems for orthogonal polynomials (see [8] for an introduction and bibliography 
of earlier works, and (23 El EBl ES EQ El ESI CES] ) • 

Note, in particular, that detY(z) = 1. (From the conditions on Y(z) it follows that 
dety(^) is analytic across the real axis, has all singularities removable, and tends to the 
identity as z — > 00. It is then identically 1 by the Liouville theorem.) The solution is unique. 
Indeed, if there is another solution Y(z), we easily obtain by the Liouville theorem that 
Y{z)Y{z)- 1 = 1. 

The case of the weight w(x) = \x\ 2a exp(— x 2 ), a > —1/2, was considered in [27] (the 
argument is straightforward to generalize to our case). Namely, [27J gives us a ready proce- 
dure to calculate asymptotics of Y(z) in this case (by applying a series of transformations 
to the problem (a)-(d)). By (fT2"|) . these results are then interpreted as asymptotics of the 
polynomials p n (z) and their Cauchy transforms. 

In the next section we shall derive an expression for (d/daj) In D n (ai, . . . , a m ) in terms 
of the matrix elements of (fl2|) . and in the section after that compute the asymptotics of 
Y{z). 
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The existence of the system of orthogonal polynomials Pk(z) = XkZ k + ■ ■ ■ satisfying (jSj) 
with nonzero leading coefficients x& for atj real, atj > —1/2, is a classical fact. Moreover, the 
coefficients k\ = Dk/Dk+i are regular functions of all ctj (as follows from the determinantal 
representation for Dj- [32], Ej). For all complex aij in any fixed closed bounded set of the 
half-plane Sfax,- > —1/2, j = 1, . . . ,m (denote the set of such m-tuples {ozj}jL 1 by Cl), we 
shall prove below the existence of a solution to the Riemann-Hilbert problem for all n larger 
than some no > 0. Its asymptotics will be explicitly constructed. We shall also see that the 
coefficients x&, k > no are nonzero and finite for all such ocj. For k < no, the coefficients 
k\ are regular and nonzero (as follows from the determinantal representation) and are 
nonzero outside of a possible subset Cl of Cl. When all aj except for <x, (for some jo) are 
fixed provided only {a J }™ =1 G Cl, we denote the set of values of <x, such that {a,,}™^ G Cl 
by Cl(ai, . . . , (Xj -i, *, &j +i, • • • , cx m ). As a consequence of the determinantal representation, 
the system of the orthogonal polynomials exists, and the formula ([9]) holds for aij in the 
jo-component of Q, outside the set £l(ai, . . . , aij -i, *, <Xj +i, . . . , a m ). Throughout Section 3 
(and, hence, in the differential identity obtained there) we assume that {aj}™ =1 ^ Cl. (There 
is no such condition on aj in Section 4.) This provides in Section 5 a proof of Theorem 1 for 
{«j}j=i outside the set Cl. However, as we shall see below, the error term in the asymptotics 
of D n is uniform for all {aj}^=i G Q. Theorem 1 will follow then by continuity. 



3 Differential identity 

Throughout this section, we consider n a fixed positive integer and {a.,}™^ G Cl \ Cl (see 
Section 2). Let us fix some v from 1 to m, and denote a u = a. The orthogonality property 
of the polynomials p n (x) = x n x n + • • • implies that 

S 

Pk{x)x : 'w{x)dx = —, j = 0, 1, . . . , k, k = 0,1,2,... (13) 



Using ([9]) and this relation, we have 

, , n— 1 Ti—l / n—1 

3=0 j=0 ' <3 j=0 



— lnD n (ai,...,a m ) = — lnj^xj 2 = -2^-^ = -2^ / p j (x)p' ja (x)w(x)dx 
a a 3=0 j=o Kj j=o 

/oo /n-1 \' 
| ^JPy(^) J w{x)dx. 
00 Vi=o / a 

(14) 

Here the prime and the lower index a stand for the derivative w.r.t. a. Below we also use 
derivatives w.r.t. x denoting them with the prime and the lower index x. 

Note that p n (x) are analytic functions of aj as follows, e.g., from their representation as 
a determinant and regularity of the leading coefficients for {otj}f =l G Cl \ Cl. 
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By the well-known Christoffel-Darboux formula (e.g., [32]), 

n-1 

= ^-lfpli^K-iW - p n {x)p' n _ lx {x)), (15) 

3=0 

where bj are coefficients in the recurrence relation for our polynomials: 

bj-iPj-i(x) + (cij - x)pj(x) + bjp j+1 (x) = 0, j = 1, 2, . . . 
Let us fix the notation for the 3 leading coefficients of the polynomials p n (z) as follows: 

Pn (z)K- 1 = Z n + f3 n Z n ~ l + ln Z n ~ 2 + ■■■ 

Comparing the coefficients at the powers z n+1 , z n , and z n ~ l in the recurrence relation, we 
obtain the following identities we shall use later on: 

bj = ^- t a j = /3 j -p j+1 , (%=±) =7i-7i+i-^ + ^i+i- (16) 

Substituting i jTST) into pij) . differentiating it w.r.t. a, using the orthogonality and the 
above expression for 6 n _i, we obtain: 

— ln£) n (ai,...,a m ) = -n 1 (Ji - J 2 ), 

da x n _i x n , , 

/oo ;>oo V ' / 

i^K-LxW^W^, ^2 = / p' ntX {x)p' n _ liCt (x)w(x)dx. 
•oo — oo 

Note that since (IT7|) contains polynomials with indices n — 1 and n only, it is already 
an identity one could use to obtain an asymptotic expression for the logarithmic derivative. 
This identity is valid for a general weight w(x). However, to obtain an asymptotic expression 
for (d/da) ln.D n , one would have to integrate asymptotics of the polynomials over the whole 
real line, which would make the calculation rather cumbersome. We shall now see that the 
particular structure of our weight w(x) allows us to simplify fJTTjl considerably and reduce 
it to an identity involving only some polynomial coefficients and the values of polynomials 
and their Cauchy transforms at the points fij. 

Let us evaluate J\. Choose points Ck, dk according to the inequalities: 

- oo = ci < /ii < di = c 2 < fi2 < d 2 = c 3 < • • • < d m -i = Cm < fi m < d m = oo (18) 

For a small e > 0, split the integral into 3m ones over the intervals (c^, e), (/i^— £, /ifc+e), 
(fik + e, dk). The middle ones can be written in the form 

\ x ~ Hk\ 2oik f{x)dx, 



where f(x) is a smooth function on — e, Hence, this integral is of order 0(e 2ak+1 ), 

small for 3?a fc > —1/2. For the same reason 

f(x)\x - ^\ 2a % k k -e = (ffrk + e)- f(jM k - e))e 2a > = 0( e *» +1 ). 
Now we can integrate each of the remaining integrals 

J Uh+e 



by parts and then take the limit e — > to obtain J x . Since w(x) decreases ase 1 at infinity 
and 

\j=i x Pi J 

integration by parts gives (we again use the orthogonality to write ^ Vn-\P" n x a wdx = 



■J m r c 

j 1 = _ n ^L_2Va J / 



X 1 ."\ roc p/ />oo 

x - f p n _i(x)^^ w(x)dx + 2 / p n _i(x)p^ Q (x)xit;(x)£Zx, (20) 

> X — fij J-oo 



where the first m integrals are taken in the principal value sense. Note that they converge 
for "Raj > —1/2. This can be seen directly as follows. Take the fs one and write it in the 
form 



' X ^ f(x)dx = lim ( r -(jij - x) 2a ^ 1 f(x)dx + [ (x- pL j ) 2a ^ 1 f(x)dx) , 



oo x A*j \J-oo Juj+e 



where f(x) is a smooth function at x — fij. Change the variables y — fx 3 ■,— x in the first, and 
y = x — fij in the second integral on the r.h.s. Then the sum reduces to 



/•oo 

J y 2 ^\f{ H + y)-f{ N -y))dy. 



The integrand is of order 0(y 2aj ) for small y, and therefore the integral converges as e — > 0. 
We now simplify the integral. Adding and subtracting p' na (fij) gives: 



/ ^Pn,a\ X ) / \ i 

Pn-il^J u>(£)aa; = 



/ 

J -oo A*J 
/■°° ( ,Pn,a( X ) - Pn,a(t*j) / \ 7 , / / \ [°° p n -i(x)w(x) 

/ Pn-iW w(x)da: + p n)a (^) / eta. 

The ratio in the first integral on the r.h.s. is obviously a polynomial of degree n — 1 in x — fij 
with the leading coefficient x/ na . Therefore, by the orthogonality, 



/ \Pn,a( x ) Pn.ai^j) , \ , ^n.a 
Pn-l{ x ) W{X)dx = . 

8 



The last integral in the expression (120!) for J\ can be rewritten as follows: 

p n _i(x)(x' a x n+1 + (x n (3 n )' a x n + (x„7 n )' a x ri ~ 1 H )w{x)dx. 



Expanding here x n+1 and x n in terms of the polynomials Pj(x), we obtain by the orthogonality 
and (fl6l) that this integral equals 



x„ 



Xn-l 

We, therefore, finally obtain 



xi 



2€n \ %i 



Xn-l 



(21) 



Ji = -(n + 2 V aj)-^- - 2 V ajp' (fij) \ 



p n -i(x)w(x) 



x — 



dx+ 



x„ 



■J / \ 2 

X„ „ / X n __i 



(22) 



The analysis of J2 is similar (and simpler). We obtain that 

p n (z)w(x) 



x 7 - m 

J 2 = 2 n ~ 1,a - 2 OijPn-l^j) 



x„ 



00 A^i 



(23) 



again with the principal value integral. 

Note that at aij = in the above two equations the terms multiples of cej disappear. 

Thus, by (TT7]). f[2"2"j) . f[2"3"j) . and ffT2l . we can now write our identity in terms of the matrix 
elements of Y(z). Some care is needed when comparing the Cauchy transforms in the second 
column of Y{z) with the principal value integrals in J\ and J2: in particular, note that Yu(z) 
and Yti{z) are unbounded at z = fij if Jtaj < 0. 

Consider —1/2 < dtatj < 0, atj 7^ 0. The contribution to the matrix elements Y^z) and 
Yti{z) containing a singular part has the form 



dj \x-pij\ 2a i 

x — z 



f(x)dx 



(24) 



with a function f(x) analytic in a neighborhood of x = fij (analytic continuation off the real 
axis). On the other hand, the corresponding term in (122]) or (123]) has the form 



|2q,- 



-f(x)dx 



(25) 
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1 




• 

z 








te- 



c 



Figure 1: Integration contour for (J2 



with the same f(x) and the integral taken in the sense of the principal value. The difference 
of (j24p and (j251) can be written as the principal value integral 



[z - Hj) lim 



t ' ?f{x)dx + / . 1 r/(s)ds ) . (26) 

[x -z){x- fij) J N+e [x -z){x- fij) 



Consider 



(x — fij) 2ai 



-f(x)dx, 



>c (x-z)(x-Hj)' 
where the contour C is shown in Figure 1. 

On the one hand, this integral equals its residue at x = z 

2m(z-^- 1 f(z). 



(27) 



{21 



On the other hand, we can rewrite the integral as a sum of parts along the real axes, a circle 
(of radius e) around x = fij, and the rest of the contour (denote this rest Co): 



(x — fij) 2aj 



f(x)dx 



C (X- Z)(X- fij) 

£ 2oy e -2ma>j _ g 



{e 2ni aj _ e -2nia j) j \ X W _ f ^ dx+ 



(x — z){x — n j) ' 



2a j p -2maj _ p 2niaj Q/ £ 2 aj +l^ 



(29) 



Hj-z 



2aj 



+ 



z-fij 



+ Fi(z), 



where F\(z) is regular as z — > fij and given by the integral over Cq. The last 2 equations 
give the first integral in (1261) : 



(x-z)(x-Hj) 



f(x)dx 



71 



iz-^y^m- 



£ 2aj 



F 1 (z)+0(s 2 ^+ 1 )/(z-fi j ) 



sin27raj V rj ' J v ' 2a j(z — fij) J v 2isin27r<x,- 
The second is obtained by a similar analysis with the cut of the root now stretching to the 
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right hand side of \ij and x — fij in the numerator replaced by (x — Hj)e 171 '. We obtain: 

\2ai 



\X — 



(x — z){x — jlj) 



f(x)dx 



sin27raj 2aj(z — fij) 2isin27ra 



where -^2(2) is regular 

Now ( 1261) can be written as follows (note, firstly, that e 2a3 -contributions from the small 
circles around fij from the 2 integrals cancel each other, and secondly, that F\ 2(2) (z — fij) = 
0{z-n s )): 

d 3 \ r — i, .|2Q!j I'd] I™ _ .. \2a-j 

1 ^ f{x)dx- / 1 ^ f(x)dx = 



x z Jc 3 X ^3Q^J 

cosnaj 

as z — > fij, Qz > 0, —1/2 < 9ftaj < 0, ctj 7^ 0. Thus, the difference of Yj 2 (z), j = 1,2 and 
the principal value integrals has a nonvanishing part at fij given by the first term in the 
r.h.s. of this equation. The difference vanishes for tftotj > (as a similar analysis shows, this 
difference has only (z — fj,j) m ln(z — /ij)-type singular terms with m > 2k + 1 for <x,- = k + 1/2, 
= 0, 1, . . ., z — > fij, and no irregular terms for other positive values of aj). Finally, there 
is a constant difference for aj = 0. 

We denote by Y vp ([ij) the matrix Y with the integrals of the second column replaced 
by the principal value integrals. Namely, 

Y^(^)= \im(Y k2 (z)-S k (z)), (31) 

where the limit is taken along a non-tangential to the real line path in C + , and 

f e ~ i7Taj ^( z ~ I'/f'-hU',)- for -1/2 < ftaj < 0, Uj + 
Sk(z) = I nifkM, ' for ay = (32) 

{ 0, for fcaij > 

with 

fi{Hj) = {2Txix n Y l p n (^ j ) Y[ \fJij - /i s | 2 ° s exp(-//|), 

fl{V>i) = -Xn-lPn-l(^j) JJ \Hj - Us?** exp(-^ 2 ). 

s+3 
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We then have by (PH), ([22]), ([23]), and (fT2|), 
^lnD n (ai, . . . ,a m ) = -n(mx n _i)' a - |n + 2^a,^ (lnx n )' Q - 2 (^^j ( ln ~^") + 

2 £ a, (x-^x^n^-));^?^) - x n _ 1 (x-\F 21 ( / i,))' a r i v 2 p (/i J )) + 2 [7^ - ■ 
3=1 

(33) 

This can be somewhat simplified by writing out the derivatives of x n Yn(/Xj), x~ x 1 Y2i(n 3 ), 
and using the identity 1 = det Y(ji) = det Y vp (/j,) = K u (/i)K, v 2 p (/i) - F 21 ( / u)F 1 v 2 p ( / u). 

We obtain then, recalling that a = a u , the following 

Proposition. Let aj, dlotj > —1/2, j = 1, ...,m be such that the system of orthogonal 
polynomials Pk(z) with finite nonzero leading coefficients satisfying exists. Fix n > 1. 
Let p n = x ra (x n + (5 n x n ~ x + 7„x n ~ 2 + ■•■), the matrix Y be given by KW\) , and Y vp be defined 
in (E2P- Then 

ln-D w (ai, . . . ,a m ) = - + 2^07^ (In x n x ra _i)' Qi/ - 2 (^-77^) + 



d 



da u 

m 

2$>, {YjiM'cJ&M - Y-iWaY-TQij) + (lnx^^o^yu^on?^)) + 

3=1 

2 b'n,a„ ~ PnP' n ,a J , V = 1, . . . ,m. 

(34) 

Note that x^_ x = lim^oo ^i-i an d 7 n , /3 n are expressed similarly. Thus we can find 
the r.h.s. of f[34l) asymptotically as n — > oo provided the asymptotics of Y(z) are available. 
These are found in the next section. 



4 Asymptotic analysis of the Riemann-Hilbert prob- 
lem 

In this section we are guided by [27[ [13 where the necessary steepest-descent analysis of the 
problem (a)-(d) of Section 2 was carried out. We are now faced only with a straightforward, 
although cumbersome, calculation. We shall see that in order to obtain asymptotics of 
the r.h.s. of ( !34l) to the needed accuracy, we have to compute the first two terms in the 
asymptotics of the coefficients x n , /3 n , 7„; and only the main term in the asymptotics of 
Y(nj). In this section we assume only the condition Statj > —1/2 for 07. 

For the rest of the paper, we assume that 

A i G(-l,l), j = l,...,m, Xj^Xk if j ^ k - 
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4.1 U, T, and S transformations of the Riemann-Hilbert problem 

As usual, we perform a series of transformations of the initial problem for Y(z). The first 
one Y — > U is a scaling: 



(35) 



The second one C/ — > T is given by the formula 

U(z) = (2n)- 4<T3/2 e nZ,T3/2 T(2)e n(9{z) - i/2)<T3 (2n)-- 4,Ta/2 , (36) 

where 



A = ^2a j: Z = — 1 — 2bi2, 



/! 2 
ln(z - s)^(s)ds, zeC\(-oo,l], ^(z) = _Vl - z 2 . 
1 7T 



(37) 



Below we always take the principal branch of the logarithm and roots. The function ip(z) 
is the scaled asymptotic density of zeros of p n (z). (Note that J^_ 1 ip(x)dx = 1.) The factor 
e ng(z) can th ere f ore De regarded as a rough approximation for the polynomials. The function 
g(z) has the following useful properties: 

g + (x) +g-{x) - 2x 2 - I = 0, for x e (-1,1) 

g+(x) +g-(x) - 2x 2 - I < 0, for iel\[-l,l] 

( 2m, for x < -1 ( 38 ) 

0+0*0 - #-0*0 = { £ for a: G [-1, 1] 

I 0, for x > 1 

From the Riemann-Hilbert problem for Y(z) we obtain the following problem for T(z): 

(a) T(z) is analytic for z G C \ R. 

(b) The boundary values of T(z) are related by the jump condition 

T + (*)=T4*)[ Q A ^S,(.)-p-(i)) )' -e(-l,l)\U- 1 {A,}. 

/l n m It - A .\^ ip n{g+{x)+g-{x)-2x 2 -l) \ 

T + (x)=T_(x)^ Q ll i= llX ^ \ ), iGR\[-l,l]. 

(39) 

(c) T(z) = I + 0(l/z) as z -> oo. 
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A 




Figure 2: Contour for the Riemann-Hilbert problem (m = 2). 

(d) Behavior of T(z) as z — > Xj is the same as for Y(\/2nz). 

Note that this problem is normalized to 1 at infinity, and the jump matrix on (— oo, — 1) U 
(1, oo) is exponentially close to the identity (see (138}) ). We have to exclude small neighbor- 
hoods of the points —1 and 1, where g + (x) + g_(x) — 2x 2 — I is close to zero, to have a 
uniform bound. 

Now let h(z) be the analytic continuation of 

h(x) = g + (x)-g-(x) = 2vre 3 ^ 2 J%{y)dy (40) 

to C \ ((— oo, — 1] U [1, oo)). Then a simple analysis shows that $lh(z) > for Qz > 0, and 
dth(z) < for 'Sz < in some neighborhood of (—1,1). We again exclude neighborhoods of 
the points —1 and 1 for a uniform estimate. 

It turns out it is possible to make use of this fact and split the contour on (—1, 1) as shown 
in Figure 2., transforming again the Riemann-Hilbert problem accordingly (steepest descent 
method of Deift and Zhou). Then more elements of the jump matrix become exponentially 
small in n, and hence asymptotically negligible, and what remains can be reduced to matrices 
with constant (in z) elements. These jump matrices are different in the m + 3 neighborhoods 
U±i, U\. (discs around the points ±1 and Xj of some sufficiently small fixed radius 5), and 
Uoo = C \ (Ui U U-x WJL. 1 U\ } ). Note that in U\ j we will have to add an additional cut 
roughly perpendicular to the real axis (because of nonanalyticity of the absolute value). Our 
task is then to construct approximate solutions (parametrices) P±i(z), P\.(z), Poo(z) in the 
corresponding neighborhoods and match them on the boundaries dU±i, dU\.. The solution 
in Uoo is given in terms of elementary functions; in U±±, in terms of Airy functions [12J; 
and in U\., in terms of Bessel functions [27J [31]. The "matching" is performed using the 
known asymptotics of these functions at large arguments and one more "i?" transformation 
of the Riemann-Hilbert problem. This allows us to compute any number of terms in the 
asymptotics of Y(z). 

We start by splitting the contour (—1,1) into lenses (Figure 2). Define an analytic 
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continuation of YYjLi \ x ~ A 4 ?! 2 " 3 as follows: 

f nr=i(A, -^, for 9te < Ai, 

^) = < n;=i(^ - a,) 2 ^ nr=, + i(A, - *) 2ay , f ° r ^ < ^ < w ^ = 1, . . . ,m - 1 

lnr = i(^-A,) 2Q % for3te> A m . 

(41) 

Splitting of the contour is possible because of the following factorization property of the 
jump matrix on (—1,1): 

e -„h(x) uj(x)\ = [ 1 OW 1 0\ 

e"'* 1 ' J ^(x)-^"^) lyy-^x)" 1 y ^(rc)- 1 e- nfc ( x ) l/' 1 ' 

Define a new transformation of our matrix-valued function as follows: 
' T(z), for z outside the lenses, 

T(z) ^ ^^-lg-nh^z) ^ ' for z in the upper part of the lenses, 



S(z) = < 



(43) 



/ ( : > ( ^^-i e nh(z) i J -i for z in the lower part of the lenses. 



(Note that we had to contract the lenses at Xj because we want S(z) to be analytic in each 
part of the lenses.) 

Then the Riemann-Hilbert problem for S(z) is the following: 

(a) S(z) is analytic for z G C \ E, where S = MU U*™+%. 

(b) The boundary values of S(z) are related by the jump condition 

S+(x) = S-(x) ^u^-lgfnhix) -^j > X e U j=l" 2 ^j) 

where the plus sign in the exponent is on £ 2 j, and minus, on £2,7-1, j = 1, 2, . . . , m + 1, 

M^=S-(*)f ^ "j^), xe(-l,l)\u- 1 {A J }. 



S+(x) = S_(x] 



(I {) 1 J, xeR\[-l,l]. 



(44) 



(c) S(z) = I + 0(1/ z) as z -> 00. 

(d) For Sfto,- < 0, the matrix function S^z) has the following behavior as z — > A^: 

S(z)=o(\ \lZxP>)> asz->A,-,zeC\E. (45) 
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For dtaj > 0, the matrix function S(z) has the following behavior as z — > Xj: 

O ( 1 J , as z — > Aj from outside the lenses, 

^ ■--) = <; /u_\.ii4 n ( 46 ) 

O ( ^ J j_2a^ -j) > as z — > Aj from inside the lenses. 

Recalling the remarks above, we see that, outside the neighborhoods U\., U±\, the jump 
matrix on E 3 -, j = 1, . . . , 2m + 2 is uniformly exponentially close to the identity. 

Let us now start constructing parametrices which give (in their respective regions) the 
leading contribution to the asymptotics. 



4.2 Parametrix in U, 



3G 



First, we expect the following problem for the parametrix in U^: 

(a) Poo(z) is analytic for z G C \ [—1, 1], 

(b) with the jump condition on (—1, 1) 



1?')' "("^^N^V. («> 

(c) and the following behavior at infinity 

F^z) = J + oQj , asz^oo. (48) 

A solution P OQ (z) can be found in the same way as in |28j . 

i / -i ■/ _1\ \ / 1 \ 1/4 

1 Nff, / a + a 1 — z a — a 1 ) \ s _ m , N / z — 1 x 



P «W = 5 (P ~ ) "(i(a- a -') a + a"> J^W"" °« = (iTT J ' (49) 
where the cut of the root is the interval (—1, 1). (Note that det Poo(z) — 1.) Here 

v 7 ^^! f 1 lnu;(0 d£ 



= exp 



2tt yr^e^-e 



Poo = lim X>(z). (50) 



The Szego function T>(z) is analytic outside the interval [—1, 1] with the boundary values 
satisfying V + {x)V_{x) = lu(x), x G (—1, 1) \ U™ =1 {\j}. We can calculate the integral in (j501) 
directly as follows. First, replace co>(x) by u;(z/x); second, take the derivative of the integral 
w.r.t. v. Set then £ = cos and integrate in the complex plane of e ld . After that integrate the 
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result w.r.t. v over (0, 1), and finally, use the value of the integral at v — 0. This procedure 
is actually a simple "scalar" analogue of the whole present work. For \j G (—1, 1), we obtain 

m m 

V(z) = (z + Vz^~l)- A l[(z - \j) aj , A = Y, ( 51 ) 

Therefore, 

©oo = 2~ A . (52) 



4.3 Parametrices at z = Xj 

Let us now construct the parametrices in U\., j — 1, . . . , m. The construction is the same 
for every j. We look for an analytic matrix- valued function P\ 3 {z) in a neighborhood of 
U\. which satisfies the same jump conditions as S{z) on E n has the same behavior as 
z — > \j, and (instead of a condition at infinity) satisfies the matching condition 

P Xj (z)P^(z)=I + 0(l/n) (53) 

uniformly on the boundary dU\ j as n — > oo. 

Using the analytic continuation of ^(y), define: 

,AM = J M*)/2 = e^/ 2 n J^(y)dy, for 3* > 0, 

n J \ e - iw /i(z)/2 = e^Tr fii/>(y)dy, for 3,2 < ' 1 1 

It is easy to verify that the function is analytic outside [—1,1]. Note that h(z) was 
defined analytic outside R\ (— 1, 1). Set furthermore, 



f(z) = vr / 4>(y)dy. (55) 

Let us now choose the exact form of the cuts £ in U\. so that their images under the mapping 
C = nf(z) are direct lines (Figure 3). Note that ((z) = nf(z) is analytic and one-to-one in 
the neighborhood of C/ A , and it takes the real axis to the real axis. We have: 

C = nf(z) = n2yJl-Xj(z - A,)(l + 0(z - A,)), z -> A r (56) 

Set 

Wj(z) = f](z - \r ft (Ai - ^ x ( - if 71/2 < 1 arg/WI < 7 • (57) 

This function has the following jumps on r 3 , T 7 : 

w J+ (z) = Wj4z)e>™*, (er 3 u r 7 . (58) 
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Moreover, the functions Wj(z) and u(z) are related in U\ j as follows: 

Wj(z) = u{z)e- 2 ^ (59) 

in the region K( > 0, Ssz > n $C > 0, and in the region 3?C < 0, $sz < n $C < 0; whereas 

Wf (z) =u(z)e 2ma > (60) 

in the region 3?C < 0, Sz > n $C > 0, and in the region 3?C > 0, %z < n 3( < 0. Note 
that by definition of u(z), argcj(z) = on $tz = 0, and cu(z) is continuous through this line. 
Hence the function Wj(z) also has jumps on Ti, T 5 given by the formulae: 

H/, + (,)=H/,_ W x| e2 _^ for J Gr ;. (61) 

We look for P Xj (z) in the form: 

P x .(z) = E n (z)P^(z)W ] (z)-^e- n ^, (62) 

where E n (z) is analytic and invertible in the neighborhood of U Xj , and therefore does not 
affect the jump and analyticity conditions. It is needed and so chosen that the matching 
condition be satisfied. 

It is easy to verify that P^(z) satisfies jump conditions with constant jump matrices. 
Because of Wj(z), P^(z) has also an additional jump condition along the line ^tf(z) = in 
U Xj . Set 

P«(«) = *« j (C) = * aj (n/(2)) ) (63) 

where the cuts for ty aj (() are shown in Figure 3, and the Riemann Hilbert problem for it is 
as follows: 

(a) is analytic for (gC \ U^ =1 iy 

(b) ty aj satisfies the following jump conditions: 

*a„ + (C) = *a„-(C) (_ j J), forCGlMJTs, (64) 

*« JI+ (C) = *«,,-(C) ( e JL, J), forCer 2 ur 6 , (65) 
* ajl+ (C) = ^, . iO' ;: """ ! . forCer 3 ur 7 , (66) 
*a„+(0 = * aj ,-(0 J), forCGr 4 ur 8 . (67) 
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Figure 3: The auxiliary contour for the parametrix at Xj. 
(c) For 3taj < the matrix function ^> a . (() has the following behavior as ( — > 0: 



C-o. 



For JJofj > the matrix function (£) has the following behavior as ( — > 0: 

<MC) = 



' j^P j^j-I- ) » as z -> Aj- with C G II, III, VI, VII, 



O I ) , as Aj with C G I, IV, V, VIII. 



(68) 



(69) 



The solution of this Riemann-Hilbert problem was constructed in [31] in terms of Bessel 
functions. For our purposes, we need its explicit form only in the region II (see Figure 3). 
There, we have (see [34|): 



^ (ft - JitC \ i ^+ 1 /2iCe )e 

" V M -^-i/ 2 (Ce- m/2 )e-^/ 2 



n K aj+1/2 {(e-™ / 2 )e^/ 2 



1 

--K, 



_- Aoj -i/2(Ce" 



iTr/2\ e iTrctj/2 j ' 



CG/J, (70) 



where /,g(x) and Kp(x) are modified Bessel functions (see |T|). The solution in other regions 
can be reproduced by applying the jump conditions to (1701) . 

Taking also E n (z) from [27J,[3l], and substituting all into (1621) . we obtain for z G z(II), 
where z(//) is the image of // under the mapping ( — > z, 



1 % 



l \^ a3 {nf{z))W J {z) 



-a-3 e -n<t>(z)a 3 



(71) 
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The argument of Bessel functions is uniformly large on dU\ ■ Substituting asymptotics 
of Bessel functions for large arguments 



T />x 1/1 — i 



4( \ -2iaj 2aj 



e (-7r/4-Q!j7r/2-C)io-3 £ ^ jj 



(72) 

(and similar ones for other regions (see [271 El])) into ( 17TI) (and its analogues for other 
regions), we can verify that the matching condition is satisfied. For that we need only 
the main asymptotic term plus the error term. The fit is ensured by the choice of E n (z). 
The exponential factor in the asymptotics of Bessel functions cancels with ^- n ^ z ) a ^ in (T7T1) 
leaving a constant in z factor. (Note that f(z) = i<fi{z) — i<fi + (\j).) Moreover, considering 
further terms, we can extend (1531) into full asymptotic series in inverse powers of n. For our 
calculations we need to know explicitly the first correction term: 

P Xj (z)P^(z)=I + A 1 (z)+0(l/n 2 ) 1 

i ( a 2 ictj 



y J V 1 jV ; 2C Viaj -a) ) (73) 

W .^-03 e -( rl <l > +( X 3)+ a j'Ki/2— tu/4)<t 3 p- 1 ^), 

z G dz(II), 

where dz(II) is the part of dU\. whose C-image is in II. As the calculation for the other 
regions shows, this expression for Ai(z) extends by analytic continuation to the whole bound- 
ary dU\. (cf. [2"7l I3"l] ). Moreover, it gives a meromorphic function in a neighborhood of U\. 
with a simple pole at z = Xj. The error term 0(l/n 2 ) in (1731) is uniform on dU\.. 

Note that the absence of jumps for det ^(C) and the fact that dtotj > —1/2 implies that 
the only possible singularity of det ^(C) (at C = 0) is removable. Thus, det ^(C) is analytic. 
Moreover, the asymptotics of Bessel functions give that det\l/(C) — > 1 as C — > oo, which 
implies that det\I / (C) = 1. Using this, we easily deduce from flTTj) that also det Px^z) = 1. 

Note that it follows from the asymptotics of Bessel functions and the fact that 

D{z)W J {z)- 1 = (z + ^~l)- A e ±i ^^' ak (74) 

(where f = j or j' — j + 1 and the sign is different in different quadrants of the (-plane) , that 
all the terms (including the error term) are uniform for a& in a bounded set of the half-plane 
3?«fc > —1/2. This observation (and similar ones for dU±i) will be very important below. 



4.4 Parametrices at z = ±1 

Now let us construct parametrices in the remaining regions U±i. These are obtained by a 
slight generalization of the results of [12] (which can be viewed as the case u(z) = 1). We 
are looking for an analytic matrix- valued function in U\ which has the same jump relation 
as S(z) there and satisfies the matching condition on the boundary: 

P l {z)P^{z)=I + 0{l/n). (75) 
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The solution is: 

Pi = E(z)Q(()e- n ^ z)a3 u(z)-' 73/2 , E{z) = P^z)^)* 3 ' 2 ^™^^ (\ ) ^V 7 ™/ 12 , 

(76) 

and Q(C) is given by the expression (7.9) of [12] in terms of Airy functions (in the notation 
of P2J Q(C) = ^ CT (C))- In these formulas 

3 x 2/3 



((z) = i-ne^^z)) . (77) 

As is easy to verify, so defined ((z) is an analytic function in a neighborhood of z = 1 (the 
cut of the square root is the interval (—1,1)), and 



C(*) = 2n 2 '\z - 1) (l + l(s - 1) + O(0* - l) 2 )) 



(7f 



The argument of Airy function on dU\ is uniformly large, so we can expand it into the 
asymptotic series and proceed the same way as for dU\ y As a result we have the matching 
condition ( 1751) extended to the full asymptotic expansion in inverse powers of n. We shall 
need below only the first 2 terms: 

P 1 (z)P^(z)=I + A 1 (z) + 0(l/n 2 ) 1 

A 1 (z) = P 00 {z)«,{z)"' 2 e<*»"± ( l i[ _ 1 1 /6 ) e-™^z)-^p-\z) (79) 
z G dU x . 

The function Ai(z) is meromorphic in the neighborhood of U\ with a second order pole 
at z = 1. 

The argument for the parametrix in is similar. We just mention the solution: 

P_! = £Cz)<7 3 Q(e- <w C)a 3 e-^ CT ^(^)- CT3 / 2 , 
= P 00 (z)^( 2 ) CT3 / 2 e^/ 4 v ^ ( _^ ^ 

Here 



g-W^\(T 3 /4 e — TTt/12 1 



_ _ m /3 ~ A 2/3 . _ / h(z)/2 - ztt = e 3 -/ 2 vr £ ^(y)dy, for 3z > 0, 
W-e ^ 2 n Wj , <?,lZj ~\ e -/ i (z)/2 + OT = e 5i7r / 2 7r/f 1 ^(l/)^, for Sz < ' 

(81) 

For z — > — 1: 



C(z) = 2n^(l +z)(l- 1(1 + *) + 0((1 + ,) 2 )) . 



(82) 
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Figure 4: Contour for the i?-Riemann-Hilbert problem (m = 2). 



The first 2 terms in the matching condition are as follows: 



P^{z)P-\z) = I + A 1 (z) + 0(l/n 2 ), 
A 1 (z) = Poo{z)uj{z)^ 2 e^^ ( X ( 6 



1/6 



1 



) 



e-™^u(z)-^ 2 P^\z)^e-^()- 3 / 2 , (83) 



z g dU- 



As in ( 1731) . the error term in ( 1751) (resp., ( 1831) ) is uniform for all z G <9f/i (resp., 2; G dU-i) 
and Oj's in a bounded set. 

4.5 Final transformation of the problem 

Now the construction of the parametrices is complete, and we are ready for the final trans- 
formation of the Riemann-Hilbert problem. Let 



It is easily seen that this function has jumps only on dU±\, dU\., and parts of Ej, R\ [—1, 1] 
lying outside the neighborhoods U±i, U\. (we denote these parts without the end-points 
E out ). The contour is shown in Figure 4 (for m = 2). Outside of it, R(z) is analytic: Indeed, 
as follows from (pj), (jlBj), (|55|), (Pj), and the fact that det P Xj = 1, the function S(z)p- 1 (z), 
for example, has at most a singularity of order less than 1 at Xj, which implies, due to the 
absence of jumps, that S(z)Pr (z) is analytic in U\ ■ The argument for other regions is 
similar. 

Note that R(z) = I + 0(1/ z) as z -> 00. 




z G [Too \ E, 

z G E7 Ai \S, j = l,...,m, 



z G t/i \ E, 
2; G Z7_i \ E. 



(84) 
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The jumps are as follows: 
R + (x) = R-(x)P o0 (x)[q UKX)e 1 JPocix)- 1 , xeR\[-l-5,l + 5], 

/?-(•'■) = n-(.r)F\(.r) ( ,^_}^ hM °]i^ (a;)- 1 , x G S° ut , k — 1, ... ,2m + 2, 



u(x)~ 1 e^ nh ^ 1 



where the plus sign in the exponent is on S^*, and minus, on E^-n J = l;---> m + l> 
= i?_(x)P\ :( (^)-Poo(a ; ) _1 ? £ G £?E/\j \ { intersection points}, j = 1, . . . , m, 
= i?_(x)P±i(x)Poo(a;) _1 , x G <9£/±i \ { intersection points}. 

(85) 

Here 5 is the radius of Ui, U-\. 

The jump matrix on S out can be estimated uniformly in a as I + 0(exp(—en\x\)) , where 
e is a positive constant. The jump matrices on dUx h ±i admit a uniform expansion in inverse 
powers of n (the first 2 terms of which are given by f!73|) . (1791) . and fl83|) ): 



/ + AiCz) + A 2 (z) + ... + A k (z) + Oin-*- 1 ). (86) 

Every Aj is of order XjnK (The above expressions give us an explicit form of A\{z) in each 
of the neighborhoods.) 

We look for R{z) asymptotically in the form R(z) ~ Rq{z) + R\{z) + R 2 (z) + • ■ -, where 
Rj(z), j > 0, is of the same order as Aj. It can be shown following Theorems 7.8-7.10 of 
[T2] that for any k > 1 

k 

R(z) = R A Z ) + Oin^- 1 ), Rq = I, (87) 

j=0 

uniformly for all z and for aj in a bounded set of the half-plane dtaj > — 1/2, j = 1, . . . , m. 
The set can extend up to the boundary Jfrx,- = —1/2. 

Moreover, we can substitute this asymptotic expansion into ( 1851) and, collecting the 
terms of the same order, obtain: 

Ro+(x) + Ri+(x) + (Ro-(x) + R^(x) + •••)(/ + Ax(x) + •••), % G dU Xji±x . 

R +(x) = Ro-(x) R (z) = J, 

R l+ (x)-R 1 „{x) = A 1 (x), 

R 2+ (x) - R 2 -(x) = R 1 ^(x)A 1 (x) + A 2 (x), 

k 

R k +(x) - Rk-{x) = y~]Rk-j,-(x)Aj(x), k = 1,2, ... 

3=1 

(88) 

The main term in the asymptotics of the polynomials is given therefore by the parametrices 
at the appropriate points z. The expressions for Rk{z) follow from the Sokhotsky-Plemelj 
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formulas: 

If A 1 (x)dx If R 1 .(x)A 1 (x) + A 2 (x) ^ 

dU = dU\ U dU-i WjLi 9U\ r Note that the contours are traversed in the negative direction. 

For reader's convenience, we present a variant of the proof of ( IBTH891) . a combination of 
the arguments from [12] and [28]. First, we need a bound on R(z). Let 

A = J - I, (90) 

where J is the jump matrix for R on E# = E out U dU (see (1851) ). The jump condition and 
the behaviour of R(z) at infinity imply that 

R(z) = I + C(R-A), zeC\E R , (91) 

where 

C(f) = ±[ f(s) ( ' S 



2tt« s ~ z 

is the Cauchy operator on E#. Hence, 

R_(s) = I + C-(R-A), (92) 

where C_(/) = lining C(/) as z approaches a point s G E# \ {intersection points} from the 
— side of E#. It is known that C_ is a bounded operator from L 2 (E#) to L 2 (E#). Now 
defining the operator 

C A (/) = C_(/A), 

we represent (1921) in the form 

(I-C A )(K.-I) = C A {I). (93) 

Since by the estimates above A(s) = 0(l/n) and A(s) = 0(exp(— en|s|)) (on E out ) for 
n — > oo uniformly in Oj's (in a bounded set), and s G E#, the operator norm of Ca acting 
on L 2 (E#), ||Ca|| = 0(l/n), and J — Ca is invertible by a Neuman series for n sufficiently 
large. Moreover, ||CA(/)||L 2 (s fl ) — 0(l/n). Thus ( 1931) gives 

j R_( S )=/ + (/-Ca)- 1 (C a (1)), (94) 

where 

||i2-( a )-/|| ia(El0 = O(l/n). (95) 

Hence, by (Ell), 

= / + C[A + (/ - CaJ-^Ca^JA], (96) 
and a matrix norm for some £i > and dist(z, E#) > e±, 

\R(z)-I\ < |C(A)| + |C(( J R_-/)A)| <^ + C2 ||ii;_( s )-J|| i2(SH) ||A|| L2(SH) < | (97) 

24 



R=R 



+ 



R=R 




R 



R=RJ 

Figure 5: Deformed part of the contour 



uniformly in a and z for some ci, C2, C3 > 0, and n larger then some no- 

To obtain a uniform estimate for all z G C \ S^, we (for a z with dist(z,X#) < E\) 
deform the contour as shown in Figure 5. Here is the same as T, R with the dotted 
part replaced by the semicircle of radius E\. R{z) is defined as shown, and J is the analytic 
continuation of the jump matrix for R on E^. (It is easy to show that the continuation exists 
in a neighborhood of the original Neighbourhoods where 2 lines intersect are analyzed 
similarly.) R(z) satisfies the same Riemann-Hilbert problem as R{z) but on the contour 
The argument leading to fl96l) for R and holds for R and as well. Therefore (see 
Figure 5), 

\R(z) - I\ = \R(z) - I\ < c 3 /n. (98) 

Analysis of the analytic continuation J(z) shows that we can find the same C3 for all a/s in 
a bounded set, n > n , and all z G C \ T, R up to the boundary T, R . This is the estimate we 
need. 

We now proceed with the proof by induction. On E out , we define Aj = 0. Consider the 
function Ri(z) analytic outside J2 R , satisfying the jump condition Ri+(s) — Ri-(s) = Ai(s), 
s G Sr, and the condition R\(z) = 0(1/ z) at infinity. The unique solution of this Riemann- 
Hilbert problem is given by the Sokhotsky-Plemelj formula: 

R 1 (z) = C(A 1 ) = ± [ A X ( S )— . (99) 

2m JY.R s - z 

Because of the estimate Axis) = 0(l/n), uniform in a/s, s G T, R as n — > 00, we have 

R x (z) =0(l/n), n^oo, (100) 

uniform in oij's and z satisfying dist(z, E^j) > E\. Deforming the contour, we extend this 
estimate to a uniform one for all z G C \ S/j. By f l9~Tj) and fl9"§|) . 

fl(-z) - / - i?i(z) = C{R^A - Ai) = C((i2_ - I)A + A - A x ). (101) 

The (uniform in s G and a/s) estimates R-(s) — I = 0{l/n) f[9"81) . A(s) = 0(l/n), 
A(s) = 0(exp(-en|s|)) (on S out ), A(s) - A x (s) = 0(l/n 2 ) (on 017) imply that 

R(z) -I-R^z) = 0(l/n 2 ), n^oo (102) 
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uniformly in a/s, and z such that dist(<z, £#) > E\. By a contour deformation argument this 
result extends uniformly for z6C \ 

Next solving the Riemann-Hilbert problem for i? 2 (i.e., with the jump i? 2+ — R2- = 
R1-A1 + A 2 ) we get 

R 2 {z) = C(i?i_A 1 + A 2 ). 
By fllOOp and the estimates A^(s) = 0(l/n k ), we have as above for Ri. 



R 2 (z) = 0(l/n 2 ), n^oo, 

uniformly in a/s, z G C \ T,r. By ( 1981) . (1 102H . and the estimates for A, 

R{z) -I - Rr{z) - R 2 (z) = C(R_A - A x - Ri-A l - A 2 ) = 

C((i2_ - I) (A - Ax)) + C(A - Ax - A 2 ) + C((fl_ - I - #i_)Ai) = 0(n^) 

n — ► 00, 



(103) 
(104) 



(105) 



(106) 



with the same uniformity property in a/s and 2. 

The general fc'th induction step is carried out similarly. We have, using estimates from 
the first k — 1 steps, 

R k (z) = C /,', , A ; j = 0{n" k ). 
and for the error after k terms 

k I k j 

j=i \ j=i 1=1 



Cl(R--I) 



k-l 



r k-i 



fc-j 

R- - I - R i~ 
i=l 



C A- 



A 3 - I =0(n 



3=1 

k-l\ 



+ 



(107) 



The last 2 estimates are valid uniformly for a/s in a bounded set and z G C \ This 
concludes the proof of (I8TH89]) . 



4.6 Calculation of the asymptotics for the polynomials 

In the intersection of the region A (the area outside the lenses) and the neighborhoods U we 
have by (El), {43}, (JET 



T(z) = S(z) = P e (z)+R 1 (z)P e (z) + 0(l/n 2 )P e (z), z G U e n A, 6 = ±l,\ j ,oo. 

(108) 
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It turns out that to prove Theorem 1 using (IMj) . we need to evaluate only P\ (z) as 
z — > Aj (i.e., the main asymptotic term of T at Aj) and (1 1 8 [) as z — > oo (i.e., the first 2 
terms of T at infinity). 

We evaluate P\Az) as z — > Aj inside the region -2(11) (see Figure 3). -Paj(^) there is 
given by (TTTi and (j7Dj). 

Tracing back the transformations (11081) . fl36l) . and fl35|) . we obtain for the main asymp- 
totic term 

Y(zV2n) = (2n) nCT3/2 e nia3/2 (2n)^ 3/2 (/+0(l/n))P Aj (^)(2n)-- 4<T3/2 e n(9W -' /2)<T3 , z G z(/J). 

(109) 

For application in the identity (jM]) . we need to analyze equation (11091) as z approaches 
Xj along a path in z(II). Note first (recall fl5"Tj) ) that 

V{z)W j {z)- 1 = e iT& i{z + vG^l)- A = e ^-w(T/2-rj)(i + 0(z - A,-)), 

1 m (110) 

-i=(y/\j + 1 + \/ A j - 1+) = e^ /2 ^ )/2 , = arcsinAj- ctj = ^a k . v ; 

Therefore, the product of the first 3 factors in (171]) gives as z — > Xj: 

2~-4°"3 / e -(A+l/2)Tji-(2&j-A-aj)iir/2+n4>+(Xj) e (A+l/2)T j i+(2a j -A-aj)in/2-n(j> + (\ j ) 

y/2(l — A 2 ) 1 / 4 \ —e~^ A ~ 1 f 2 ^ T ' i ~^~ A ~ a ^^' ¥ ^ + ^ e ( ,A-l/2)T j i+(2a :j -A~ctj)i-K/2-n<i> + (\j) 

a+o{z-x 3 )). 

(in) 

By (15Ej) . we need the asymptotics of Bessel functions at a small argument for (170)) (see pQ 
for them). For a,- ^ 1/2 + fc, fc = 0, 1, . . ., 

<M0 = 

C 1 (a i )C^'+ 1 -^(C/2)-^e*f ( 2 ^ +1 / 2 )r( aj + 1/2) + C 2 ( ai )C^ +1 

(l + O(0), 

(112) 

where Cj(a.j) are constants whose precise expressions will not be important below. For 
(Xj = 1/2 + k, k = 0,1,..., the constants C\ and C3 remain the same, while C%{a.j) is 
replaced with C 5 (aj) + C 6 (a!j) ln£, and C±{a.j) is replaced with C^otj) + C&(aj) hi(. 

Now the part of the product of the next 5 factors in (ITTj) with the C 2aj terms omitted 
("regularized" part) gives for z — > Xj and ay 7^ 0: 



X 



1 "2 
1 I 



r(oy + 1/2) 
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This result will correspond to the "principal value" part of Y given by (l3Tj) . The omitted 
part of order £ 2a -? does not contribute to the "principal value" part of the total product for 
Y because all the correction terms are 0(() as ( — > and 9?Oj > —1/2. 

The expression (11 13ft is also valid for ctj = (since, e.g., (1321) for this case can be obtained 
by letting a, — > in the C 2Qj term with a, < 0). 

Finally, the product of the last factors in (1711) and (11091) gives by (1381) : 



Taking the product of all contributions, we have: 

lim Y vp {zV2^) = ( — ) {2n) Aas/2 {I + 0{l/n))2- AcT3 {l - A 2 )^ 1/4 x 



(114) 



2e 



V2 cos § (tj - tj) —j- sin § (t 3 - - Tj ) 
-iy/2 sin \{tj + ^= cos \{tj + 7}) 



V 



(2n) A / 2 T(a j + 1/2) 



tj = 2rm tjj(y)dy + irctj — 2n Y^T=j a » + ~ ^Tj), Tj = arcsin Xj. 



(115) 



(116) 



We now turn to the asymptotics of the coefficients x n , (3 n and 7™ of p n (z). Here we need 
the first 2 asymptotic terms in n. As usual, we compute them investigating the limit z — ■> 00 
of y(2). Namely, by (112). 



n-l 



lim 



2nz 



n—l ' 



U u (z) = z n + ^z n - 1 + ^z n - 2 + 
\/2n 2n 



;ii7) 



As z — > 00, we need to know asymptotics of Y(z) in the region A, which are given by the 
expressions (cf. ( 11091) ): 



Y{z^f2m) = (2n) na3/2 U(z), 
U(z) 



3 nZ<x 3 /2/ 9 \Aa 3 /2 



(2n) A ^ 2 (I + R 1 (z) + 0(l/n 2 ))P 0O (^)e 



n{g(z)-l/2)a 3 



(2n) 



-Ar 3 /2 



118) 



z e An Un. 



Let us compute R\{z) using (1591 . Consider first the neighborhood U\.. Substituting A 1 (x) 
given by (1731 into (1551 and calculating residues at a simple pole x = Xj, we obtain the 
contribution to R\ from the neighborhood U\ : 



Kl [Z) ~2m 



Aidx 



X; 



0{z 



-2< 



1 

2tx% 



Aidx 



X, 



0(z 



-2< 



V a3 o ■ 



— OijXj + sinti 



i{aj - cos(tj - Tj)) 



4n(l - A 2 ) \i(oij + cos(tj + Tj)) atjXj - sin^- 



^00 1 



(119) 
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where tj, Tj are defined in (I116p . 

To compute the contribution from the neighborhood Ui, we repeat the calculation now 
using Ai(z) from (J79l) . An additional complication is that we need to calculate residues in 
the pole z — 1 of second order. We obtain: 



p(i) _ 1 



Aicfe 1 



dUi 



X — Z Z 



- l + - + 0(z- 



-2\ \ ^oo 



1/8 -A 2 i/6 + iA + iA 2 
8n \i/6-iA + iA 2 -1/8 + A 2 



1 5Pg /-l i 
2 2 8 - 24nl « 1 



Here the dependence on ay comes from the expansion (see (IS"T1) ): 

Pi:) l-2v / 2^v / ^T + 4^ 2 (z-l) + 0((z-l) 3/2 ), asz->l. 



(120) 
(121) 



A similar calculation for U-\ gives 



R\ 



(-i) 



1 

2m 



In this case 



/ 


Aidx 1 


Jau- 


j x — Z z 


-1 
— i 


1 00 



1 \ T) a ' A 

l-- + 0(z~ 2 ) ^ 
z V V 8n 



-1/8 + ^l 2 i/6 + iA + iA 2 
i/6 - U + M 2 1/8 - .A 2 



w(z) 



1 + 2v / 2Uv / ^TT - 4^ 2 (^ + 1) + 0((z + 1) 3/2 ), 



as z — ► — 1. 



Summing up all the contributions (11191) . (11201) . and (11221) . we obtain: 



£>-<7 3 _ 

oo 



(122) 
(123) 



i? 1 = J RS 1) + J R!- 1) + ^^ ) . 
Substituting this into (I118P and using the expansions for z — > oo: 

2 ™ 



P(z) 



1 ^ 1 

we finally obtain from (1 1 1 T[) 



fit \ i ( h 



vi=i 



+ 0(2" 



#(z) = In 2; - — + 0(z 



(124) 



(125) 



n-l 



2 n - 1+ - A n-- A 
V5F(n- 1)! 



(126) 
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0n 



ctj sin t,- — a'jXj 



+ 




(127) 



n 



n 



A 



-E 



a: 



cos(i,- 



1 

An 



128) 



sin tj 



fe=i 




where = arcsin \j, A = Y^jLi a ji an d *j i s defined by (11161) and ( 1371) . The error terms here 
and in (1115p are uniform for aj in a bounded set provided only $ta>j > —1/2. Let us track 
the a-dependence of the error terms in more detail. These terms in f 1 1 2 6 [) — (1128 h are those in 
(11181) . at worst multiplied by a polynomial in atj (independent of n) and 2~- A coming from 
the expansion H125[) of D(z). The error term in fl 1 1 5 [) is the same as the one in (I109p . Now 
the error terms in (11091) and ( 11181) are those from ( JHTl) which we showed to have the above 
uniformity property. 

Thus we constructed a solution to the Riemann-Hilbert problem of Section 2 for n > no, 
atj in any bounded set of the half-plane tftctj > —1/2. By uniqueness, it gives the orthogonal 
polynomials via (|T2|) . On the other hand, the determinantal representation for the orthogonal 
polynomials shows that R(z) is an analytic function of <x,'s. Furthermore, Rk{z) are also 
analytic functions of a/s by construction. Thus, the error term in (I107p is both analytic and 
uniform in Oj's in a bounded set of the half-plane tftotj > —1/2. Therefore, it is differentiable 
in a/s (the derivative being of the same order in n and uniform in ay's). Hence, we easily 
conclude that the error terms in (11151) . f 1 1 2 6 [) — (1128 h have the same differentiability property. 
Alternatively, we could have deduced the differentiability of the error terms by noticing first 
that the asymptotic expansions of Bessel functions we used are differentiable in <x,'s. 



5 Proof of Theorem 1 

We now substitute the asymptotics (11261) — fl 1 28 [) . and f 1 1 1 5 j) into the differential identity (J34J). 
We assume first that all {a.,}™^ G Q \ Cl (see Section 2) and {0} G Q. Care is needed with 
estimation of x n . To obtain the asymptotics of x n from (11261) we need first to replace n with 
n + 1 and second, to replace \j with \j a/^y • Without this second step we would obtain a 

coefficient corresponding to the weight Yij \x — \j^j2(n + l)| 2aj 'e~ x2 . However, since the new 
Xj is inside (— 1, 1) for n large enough, and because it enters the asymptotics of >t n starting 
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with 0(l/n) term, this second replacement affects only tj in ()126p . We obtain 
-(n + 2^ % )(lnx n x n _ 1 )L-2 ^ In + 

j—l \ n / V n / a„ 

m 

2 [tU - J = (n + 2 53 «,) ln(n/2) + (2A* - l)n + 2a, + 0(l/n). 

j'=i 

Furthermore, 

m 

2 5Z"i ( y il(^)L y 2 V 2 P (^) - y 2l(^)a^l?(^) + (lnX n X n _ 1 ) , ^F 11 (/i j )y 2 V 2 P (/i j )) 



(129) 



>2> 



(130) 



a, ln(l - A*) - 2a„S^±M - 2 £ ay ln(n|A, - K\) + O fc) . 

Note that the trigonometric terms of (I126p - (I128I) . and ()115p cancel in these expressions. 

Legitimacy of differentiation of the error terms w.r.t. a (and uniformity of these terms 
and their derivatives) follows from that for the error terms in ()126p - ()128l) . and fl 1 1 5 [) ex- 
plained in the end of the previous section. 

The sum of (DSHD and flT3U]) yields: 

_d_ 

da 



lnD n (ai, a 2 , ■ ■ ■ , a m ) = (n + 2a u ) ln(n/2) + (2A^ - l)n + 2a u + a u \n(l - \ 2 V )- 



V 



v = 1, 2, . . . , m. 



(131) 



Since the error term here is uniform in all a u , we can integrate this identity. 

First, set a 2 — ■ ■ ■ — a m — 0, v — 1 in (11311) . This identity was established for a\ outside 
the set f2(*, 0, . . . , 0) (see Section 2). Note that the number of points in Q(*, 0, . . . , 0) is at 
most finite. Indeed, the function xl = x^(ai,0, . . . , 0) is a ratio Dk/Dk+i of two analytic 
functions of a\ for ^ta\ > —1/2. They are not identically zero because they are known to 
be positive for a\ = 0. Consider the function /(«i) = D n (ai, 0, . . . , 0) exp(— J Q ai r(n, z)dz), 
where r(n, ai) is the r.h.s. of (1131)) with v = 1, a 2 = ■ ■ ■ = a m = 0. Equation (1131)) is 
the statement that /'( a i) = outside 0, . . . , 0). Since /(cti) is continuously differen- 
tiable, /'( a i) = f° r all ai such that (ai, 0, . . . , 0) G Q. Moreover, as f(ai) = /(0) = 
_D n (0, . . . , 0) 7^ 0, the determinant D n (ai, 0, . . . , 0) is never zero. Hence the identity (11311) 
with u = 1, a 2 = ■ ■ ■ = a m = holds for all a± such that (ai, 0, . . . , 0) G Q. Integrating it 
over a\ from to some a%, we obtain: 

2 

a ti 
lnD n (ai,0, • • • ,0) = lnC(ai) + — ln(l - \\) + (ain + a\) In - + n(2A^ - l)«i + 

, x 2 2 , 132 , 

Inn x 



lnD n (0,...,0) + O 
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where 

C(a) = T(a + l/2)- 2a exp ( 2 / In T(s + l/2)ds + a 2 ) . (133) 



By freedom in the choice of Q, the expansion flT32j) is valid for any cci, 3?«i > —1/2. 

Second, a reasoning similar to the one above shows that equation (113ip with «i fixed 
and «3 = • • • = a m = holds for all a 2 , 5ia: 2 > —1/2. Integrating it over a 2 from to some 
a 2 gives 

9 

In D n (ai, a 2 , 0, . . . , 0) = In C(a 2 ) + y ln(l - \\) + (a 2 n + a\) In - + n(2A 2 - 1)« 2 
-2a 1 a 2 ln2\X l - A 2 | + lnD„(ai,0, . . . , 0) + O 



n 



"a 2 * n * ( 134 ) 

j=i j=i ~ j=i 



-2aia 2 ln 2|Ai - A 2 | + lnD n (0, . . . , 0) + O 



ln n 



n 



To obtain the second equation, we substituted hxD n (ax, 0, . . . , 0) from (11321) . Equation (11341) 
is valid for any a>i, a 2 , 9ftai,9fta 2 > —1/2. 

Proceeding in this way, we prove Theorem 1 by induction after m steps. 

The second representation (JH]) for the constant C (a) is easy to obtain using the following 
properties of the G-function [55] : 



G(z + l) = T(z)G(z), G(l) = l, 

f z z z(z + l) (135) 

/ \nT(x + l)dx = -ln27r- 1 ; + z\nT(z + 1) - In G(z + 1). 
Jo 2 2 

Another identity is useful for comparison with (j4j): 

21nG(l/2) = (1/I2)ln2-In0r + 3C'(-1). (136) 
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